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Abstract 

Using standard thermodynamics and previous results of the au- 
thor, this paper aims to discuss the conditions for phase equilibrium 
in a Lennard-Jones fluid. Possibilities of astrophysical applications of 
the results obtained here are discussed to some extent tj. 
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1 INTRODUCTION 



In this paper we shall explore the conditions for phase equilibrium in a fluid 
consisting of neutral atoms or molecules. Apart from being an interesting 
task in pure statistical physics, work on this problem has several astrophysical 
motivations. Fluids consisting of neutral atoms and/or molecules exist in the 
interiors and atmospheres of the giant planets and their icy satellites, but 
also in the diffuse molecular clouds which are present in large numbers in 
galaxies. A new field of possible applications of work discussed in this paper 
are the neutron stars. Recent work shows that the inner crust of neutron 
stars consists of a lattice of spherical nuclei imersed in a sea of free neutrons 
with a background of electrons [1]. 

According to |http: / / www.exoplanet.eu] at the end of December 2011., 
716 planets around other stars were known. Of all these planets, 83 objects 
have orbital semi-major axes between 0.9 and 1.4 astronomical units. On 
general grounds, it can be expected that this range of distances from a star 
corresponds to temperatures under which fluids consisting of neutral atoms 
and molecules can exist, which illustrates the planetological importance of 
the study of such a fluid. Water, as an example of a planetologically and 
biologically important molecule can exist as a liquid in a region of a planetary 
system in which the temperature is between 273-fT and 373-fT. The upper and 
lower radii of this zone (the so called habitable zone) depend on the absolute 
luminosity of a star. It can be shown that for a main sequence star with solar 
absolute luminosity, the inner and outer radii of the habitable zone are 0.95 
AU and 1.37 AU. An example of this calculation is available at the adress 
http: / /www. planetarybiology.com. 

Examples of interesting early work on the atmospheres of the giant planets 
are [2], [3] or [4]. Results of studies of such a fluid are also important in studies 
of cold interstellar clouds and protoplanetary disks which, because of their low 
temperature, contain neutral atoms and/or molecules. For an example of a 
recent observational study of protoplanetary disks see [5] . 

The first necessary step in attempting to model any kind of a physical 
system is to determine in some way the form of the interparticle potential in 
it. Logically, in order to have the system in equilbrium, this potential must 
be a combination of an attractive and a repulsive part. In the calculations 
to be discussed in this paper, the so called Lennard- Jones (LJ) potential will 
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be used. This potential has the following analytical form 




The symbol e denotes the depth of the potential,while a is the diameter of 
the molecular "hard core". Obviously,/im r _ >0 M ( r ) = oo. It can simply be 
shown that lim r ^ a u(r) = and that (du(r)/dr) = for r m j n = 2 1 / 6 <r. The 
depth of the potential well is u(r min ) = —e. The term with r~ 12 corresponds 
to the repulsive part of this potential. It is now known,for example [6], that 
much better results could be obtained by approximating this term with a 
function of the form Aexp(-Br) (the so called Born-Meyer term). However, 
the 12 — 6 potential has the virtue of simplicity. 

The following section contains a brief resume of the basic notions and 
main previous results, while the subsequent two parts contain the calculations 
and a discussion of their possible applications in astrophysics. 



2 THE BASIC NOTIONS 

Many systems, both natural and in laboratory, show a number of different 
phases, each of which can behave in a different way. The obvious question 
is what are the conditions under which these different phases can exist in 
equilibrium. The number of phases which can coexist in equilibrium within 
a system can be determined by the Gibbs phase rule,well known in theromo- 
dynamics. Coexisting phases are in thermal and mechanical equilibrium and 
can exchange matter [7]. In practical terms, this means that the temperature, 
pressure and the chemical potentials of the phases must be equal. As shown 
recently in [8] ,the chemical potential of the Lennard Jones fluid is given by 

1 ^ / 2?^ — 1 w e . 2n+l 

v = hd- bopio E ry r (^;— )(t)~ ( 2 ) 

where /j, ID is the chemical potential of the ideal gas, p ID is the pressure of 
the ideal gas,&o = 2ti<7 3 N^/3, Na is Avogadro's number, and T denotes the 
gamma function. The pressure and temperature of a L J fluid are related by 
the virial development of the equation of state (EOS): 

^ = Ea l (T)(nA s )'- 1 (3) 
kb-i l=1 
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All the symbols on the left side of this equation have their standard mean- 
ings, while on the right side a/ are the so called virial coefficients, A is the 
thermal wavelength and v is the inverse number density v = V/N = 1/n . 
The thermal wavelength is given by 

where h is Planck's constant divided by 2ti , k B Boltzmann's constant and 
m the particle mass. Due to increasing complexity with increasing order, the 
virial developement is most often truncated at second order terms. The first 
virial coefficient is 1, and the second coefficient for the LJ potential is given 
by (for example [7]). See also [9]. 

oo 

a 2 (T*)=5 ]T7,(l/n (2j+1)/4 (5) 

3=0 

where 

_ 2 j+i/2 2j-l 

v = -^jr T{ ^ ) (6) 

and 

T* = ^ (7) 

The values of the first few coefficients jj are: 7 = 1.733,7i = —2.564, 
72 = —0.8665, and the explicit expression for the second virial coefficient is: 

a 2 = ^^[1.733(-^r_ 
3 1 K k B T J 

2.56369(-i-) 3 / 4 - 0.8665(-i-) 5 / 4 - . . .] (8) 
k B T k B T 

As bo = 2ir<j 3 NA/3 is always positive, it can easily be shown that a 2 < for 
(^) > 0.32175. 

3 THE CALCULATIONS 

The second order approximation to the virial developement is 

p = nk B T{l + a 2 n\ 3 ) (9) 
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where the coefficient a 2 is given by eq.(7). Denoting two phases of a system 
by indexes 1 and 2, and applying to them the condition of the equality of 
temperatures and the pressure, one gets the following expression for the 
basic condition for phase equilibrium: 

mk B T + n\k B a 21 T\\ = n 2 k B T + n 2 2 k B a 22 T\ 3 (10) 

which can be transformed into the form 

T{n x - n 2 ) + T(n\a 2l \\ - n 2 2 a 22 \l) = (11) 

Introducing n\ — n 2 = x, this expression can be solved to give 

xi, 2 = n 1 + 7^|[1 ± [1 + 4a2 2 niA|(l + a^X 3 )} 1 / 2 } (12) 
3.1 Analysis 

Equation (12) can be expressed in the following form: 

x h2 = n 1 + A[l±(l+y) 1/2 } (13) 



where 



and 



2a 22 \ 2 



(14) 



y = 4a 22 niA^(l + a 2 in 1 \{) (15) 

Because of physical considerations, a; must be real, which implies that y > — 1. 
The special case of x = leads to an interesting result. By definition of x 
this means that n\ = n 2 , and by eq.(ll) it implies that \\/\ 2 = (a^/a^) 1 ^ 3 . 
Finallyeq. (4) leads to the result that m 2 /m 1 = (a 22 /a 21 ) 2 ^ 3 . 

It can easily be shown that y will be negative if the product a 21 x a 22 is 
negative and a 22 negative, which means that a 21 has to be positive. Using the 
analytical expression for a 2 we have shown that a 2 < for {-j^-f) > 0.32175. 
One of the conditions for phase equilibrium is Ti = T 2 = T which implies 
that ei < 0.32175A; B T and e 2 > 0.32175A; B T. 

In view of possible astrophysical applications, note that these two in- 
equalities combine a material parameter (e) and an astrophysical parameter 
(T). For example, for H 2 , which is present in the atmospheres of the giant 
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planets, e = 33.3/cb (for example [10]). This means that a 22 < for T < 10AK, 
which is a realistic limit in the atmosphere of Jupiter (for example). 
Assuming that y = —1, it follows from eq.(13) that 

x lt2 = n 1 + A (16) 

which further means that 

n 2 = -A. (17) 

The thermal wavelength is positive by definition, so eq.(17) is physically 
posible if a 2 2 < 0. The sign of a 2 depends on the material as well as the 
pressure and temperature. 

If y — a simple calculation shows that 

xi, 2 = n x + A[l ± 1] (18) 
which leads to two possibilities; either 

n 2 = (19) 

or 

n 2 = (20) 

The result expressed by eq.(20) describes a physically interesting situation. 
It shows that y = corresponds to the case in which the number density of 
one of the phases is zero. An astronomical example of such a situation could 
be the top layer of the atmosphere of a giant planet. 

The limiting case y— > + oo has potentially very interesting physical im- 
plications. If the temperature of the material is sufficiently high, A — > 0, 
y — > but A — > oo. This means that x\ t2 — > oo, which implies that n\ — > oo 
and n 2 is infinite of lower order than n\. Such a situation looks similar to 
what one could expect in the outer layer of a neutron star. 

A case which is also physically interesting is y— >■+ oo but T — > 0. This 
corresponds to A — > oo ,but A — > 0. Finally, it follows that in this case 
^i,2 = ni, which further implies that n 2 — > - the same result as the one 
obtained for y = 0. 

The function of the form [1 + y] 1 ^ 2 , assuming \y\ < l,can be developed in 
Taylor's series as 

[1+Vl -^ (l-2n)(n!)M^ ^ 
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Inserting eq.(21) in eq.(13) gives the following explicit form of the result for 

Zl,2 

x 12 = n 1 + A[l±Y i . 1 y n ] (22) 

L ^ (1 -2n)(n!) 2 4« y J v ; 

where y and A are given by eqs.(14) and (15). Assuming \y\ < 1 and taking 
into account terms up to and including y 4 ,eq.(22) reduces to: 

Xl2 = ni + A[l±(l + ^-^ + ^-^-)l (23) 
v 2 8 16 128 ;J v ; 

This approximation shows that the difference in number densities of two 
phases in equilibrium in a L J fluid depends on the parameters of both phases 
and (through the virial coefficients) on the parameters of the LJ potential 
in the two phases. The relative error 5 of this expression compared to the 
exact result given by eq.(13) is 5 < 10~ 3 for \y\ < 0.5. 



3.2 The final result 

The chemical potential of a LJ fluid can also be expresed as [8]: 

A* = Hid + 2piDd2 (24) 

where fiio and pio are the chemical potential and pressure of the ideal 
gas, and a 2 is the second virial coefficient. Imposing the condition of the 
equality of chemical potentials as a prerequisite for the phase equlibrium of 
two phases, it follows that 

Ai/Dl - HlD2 + 2(p/Dia 2 l - PlD2d22) = (25) 

Assume that the ideal gas behaves according to Maxwell-Boltzmann (MB) 
statistics. In that case, the pressure and chemical potential of the ideal gas 
are given by /j,id = kBT\n(n\ 3 ) and pid = nksT [7]. Inserting these two 
expressions into eq.(25), and introducing n\ — ri2 = x finally leads to 

ln[(l + -){^f] + 2[(n 2 + x)a 21 - n 2 a 22 ] = (26) 
n 2 A 2 

Inserting the definition of A from eq.(4),one gets: 

ln[(l + — ) x {^-fl 2 } + 2[n 2 (a 21 - a 22 ) + xa 21 \ = (27) 



6 



which is obviously valid only for n 2 > 0. 

This expression represents a theoretical form of the link between various 
parameters of two phases in equilibrium in a LJ fluid. It contains n,m,X,a 
and e for each phase. Three of these parameters (m,e and a ) are material 
dependent, while A contains a material parameter (the particle mass) and 
the temperature. In astronomical applications the number density depends 
on the nature of the object being studied. To render it applicable, various 
terms in this equation have to be expressed in a more "practical form", in 
the sense of replacing the implicit dependence on various material (or object) 
parameters in this expression by an explicit formulation. 



3.3 Possible cases 

The behaviour of eq.(27) depends on the functions A,y and a 2 which have 
been previously defined. 

In the case y = —1, using eqs.(16) and (17),eq.(27) takes the following 
form: 

3 TTX 

ln[l - 2a 22 \\x\ + - ln[— ] + 2(ma 21 + Aa 22 ) = (28) 
2 m\ 

under the condition 2a 22 \ 2 x < 1. 

If y = 0, there are two possibilities for x, given by eq.(18): x = ri\ + 2A 
and x — ni. In the first case, it can be shown that 

" l|1 + i x( - + if )l + 

3 , r I?l2 n r/ 1 

- In — + 2 x (m + rg )a 2 i + 

2 mi a 22 \l' 

n 2 x (a 21 - a 22 )\ = (29) 
In the case x — rii one gets the following form of eq.(27): 

ln[l + — ] + 2 x [n 2 (a 2l - a 22 ) + ma 2 i] + - ln[— ] = (30) 
n 2 2 mi 

The difference a 2 i — a 22 is material dependent, through the values of e and 
a. 

An interesting limiting case corresponds to y— >+ oo. If A and rii are finite, 
this case corresponds to rr 12 = rii — n 2 — > — oo. Physically speaking,this 
means that the density of one of the phases is arbitrarily high and tends 
to infinity, while the other phase has a finite value of the density, which is 
remindful of the situation encountered in neutron stars. 
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4 THE PRACTICAL FORM 



The term containing the logarithm can be developed into series to give: 

* ™(-l) l ,x 



ln|0 + ^)]=i:^-(^)' +1 
n 2 i—Q t + J- ^2 



(31) 



which is convergent for \(x/n 2 )\ < 1 (that is,|ni — n 2 \ < n2)and the explicit 
form of the function x is given by eq. (12). Introducing the definition of x, this 
crietrion of convergence amounts to n\/n 2 < 2. The relative error of eq.(31) 
pushed to 3 terms compared to the exact result for the logarithm is 8 < 10 -2 
for \(x/n 2 )\ < 0.5. Using eqs.(5)-(7) for the second virial coefficient, it can 
be shown that 



2 00 _2-? +1 /2 2j — 1 

a 2 i - a 22 = --kN a oI ^ — T[—^—\ 

(/3ei) ^ -(—Y((]e 2 ) 4 



(32) 



where (3 = 1/ksT and the second virial coefficient is approximately given 
by eq.(8) Limiting the development in eq.(32) to terms up to and including 
j — 2, it follows that 

a 21 - a 22 3.62959iV j 4^(/3ei) 1/4 x [1 - 

5.369387V A( 7 1 3 (/?e 1 ) 3 / 4 [l - (^) 3 (^) 3 / 4 ] - 

<7i ei 

1.81488A^x 3 (/3 ei ) 5 / 4 [l - (^) 3 (^) 5 / 4 ] - . . . (33) 

A practical problem related to the applicability of eq.(32) is the estimate of 
the relative error of this series, taking into account that there is no solution 
in closed analytical form for 021 — 022- It can be estimated by choosing 
the parameters in eq.(32),and then determining the relative difference of 
the developement with j and j + 1 terms. Taking arbitrarily that j = 3, 
a 2 = 2<7i,/3e 2 = 0.05, and that /3ei G (0.001,0.9), it follows that the relative 
error is 5 < 0.025. 
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Inserting eqs.(5)-(7),(13),(31) and (32) into eq.(27), it follows that 



" (-1)' n 1+ A[l±(l + y)V 2 ] 

2 o 

(3/2) ln(m 2 /mi) + 2 x n 2 x -irN A af 

O 



J"=0 

s 

X 



4j! L 4 



2j+i /■ '2\3/ o \2i±i 



(^)^-(-m 2 )^ 

+2x(n 1+ A[l±(l+rtx^ 

00 _2i+(V2) 

^ 47" 

j=o ^> 

xr[(2j-l)/4](/3 ei )^' +1 )/ 4 

= (34) 

This equation links various parameters of two phases in equilibrium of 
a LJ fluid, and it can be explicitely written with an arbitrary number of 
terms. Taking only the terms with j G (0, 3) and / G (0, 3) in various sums in 
eq. (34), expanding out all products and integer powers leads to the following 
result 

(3/2) \n[m 2 / mi } + — x [1 - (1/2)— + 

n 2 n 2 

(l/3)(-) 2 - (l/4)(-) 3 ] + N A n 2 at(l + ±) 
n 2 n 2 n 2 

x [7.25918(/3 ei ) 1/4 - 3.62959(/3 ei ) 5/4 - 

10.7388(/3 ei ) 3/4 - 1.7898(/3 ei ) 7/4 ] 

+iV A n 2 a 2 3 (/3e 2 ) 1 / 4 [10.7388(/3e 2 ) 1 / 2 - 

7.25918(1 - (l/2)/3e 2 )] = (35) 

where x has been defined in eq.(13). 

A simplifying case of physical interest is j3 — > 0, which corresponds to high 
temperature, when several terms in eq.(35) tend to zero. Another interesting 
simplification corresponds to x — 0. The last two expressions can be further 
simplified by assuming that m 2 = mi, which implies that also A 2 = Ai 
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5 POSSIBLE APPLICATIONS IN ASTRO- 
PHYSICS 

Equation (34) links two kinds of parameters, and it can be applied in two 
ways. It contains the number densities and the thermal wavelengths of the 
two phases, but it also implicitely contains the second virial coefficients in 
the two phases. This expression can be applied to several kinds of systems: 
cold diffuse clouds in our galaxythe interiors and the atmospheres of the 
giant planets or icy satellites, and possibly also to neutron stars. 

A cold cloud sufficiently far from any star or other source of radiation 
certainly contains neutral atoms and/or molecules, which are known to inter- 
act via the LJ potential. Immagine that such a cloud contains two phases in 
equilibrium,which means that their pressure,temperature and chemical po- 
tential are equal. Take, for example, that phase 1 is the gas at the outer edge 
of such a cloud, and that phase 2 is the gas some distance below the top. If 
the particle number density in one of the phases, ni, can be measured, eq.(34) 
provides a theoretical answer to the following question: What is the value 
of ?7-2 so that the cloud as a whole remains in phase equilibrium? A possible 
consequence of the impairment of phase equilibrium can be the onset of star 
formation. 

Globally speaking, giant planets and the icy satellites consist of fluids 
and deep in their interiors they may contain a dense central core. The fluid 
envelopes are made up mostly of hydrogen and helium, with (possible) small 
additions of heavier chemical elements. The cores are presumed to consist 
of combinations of refractory elements. The exact combination of elements 
making up a planetary core can not be determined directly. It can only 
be checked indirectly, by constructing a model of the interior structure as- 
suming some composition of the core, and then comparing the observable 
consequences of a model with observed data. 

In the outer layers of the atmospheres of the giant planets matter is neu- 
tral,and can accordingly be described by a LJ potential, while in the deeper 
regions ionization, strong ion coupling and electron degeneracy become im- 
portant [11]. Modelling the interior structure of such a system is a complex 
task in statistical physics; it demands the knowledge of the chemical composi- 
tion and the equation of state. The chemical composition can be determined 
from observation only for the uppermost layer of a giant planet or a satellite. 
Regardless of the exact form of the equation of state, pressure and density of 
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matter increase with depth. Matter making up a planet or satellite consists 
of atoms and molecules, which under increasing pressure become excited and 
finally ionised. 

Taking these considerations into account, how can eq.(34) be applied? 
Material parameters of the highest levels of the atmosphere can be de- 
termined from observation. In order to apply the calculations discussed in 
this paper,it has to be assumed that that the LJ potential is applicable at 
and/or near the upper layers of the atmosphere of the object, and that the 
object is isothermal, at least until a certain limiting depth. The next step 
is to determine, within any particular theoretical framework, the number of 
phases which exist within the isothermal region of the object being consid- 
ered, and assume that these phases are in equilibrium. When this number 
is known,eq. (34) gives the possibility to calculate the values of the second 
virial coefficient 022 and of the thermal wavelength A2 in the phase immedi- 
ately below the surface. It becomes possible by repeating this procedure to 
"follow" changes of the second virial coefficient and the thermal wavelength 
with depth below the surface. 

Interesting conclusions can be drawn concerning the time evolution of 
the equilibrium in a LJ fluid. Immagine that in a given moment of time 
one of the parameters in eq.(34) changes. A change easiest to envisage is 
the change of the particle mass in one of the phases. One would expect that 
in such a situation the equilibrium would be disturbed and that the system 
(or the part of it which was considered) would start evolving in time. There 
is however another possibility: that some other parameter (or parameters) 
changes in such a way that the equilibrium is maintained. Details will be 
discussed separately. 

6 CONCLUSIONS 

In this paper we have analyzed the phase equilibrium in a Lennard Jones 
fluid. The starting points were the basic conditions for phase equilibrium, 
well known in statistical physics, and the recent expression for the chemical 
potential of a Lennard Jones fluid. The result of the calculation is an ex- 
pression for the condition for phase equilibrium, given by eq.(27) in its most 
general form,or by eq.(34) for the case of a Lennard Jones fluid. In deriving 
eq.(34) it was assumed that the ideal gas contribution to eq.(25) is governed 
by Maxwell- Boltzmann statistics. It would be interesting to repeat the same 
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calculation, but with another form of statistics. Several possible applications 
in astrophysics were indicated, including one unexpected - the neutron stars. 
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